eDNA based bycatch assessment in pelagic fish catches.

Pelagic fish like herring, sardines, and mackerel constitute an essential and nutritious human food source globally. Their sustainable harvest is promoted by the application of precise, accurate, and cost-effective methods for estimating bycatch. Here, we experimentally test the new concept of using eDNA for quantitative bycatch assessment on the illustrative example of the Baltic Sea sprat fisheries with herring bycatch. We investigate the full pipeline from sampling of production water on vessels and in processing factories to the estimation of species weight fractions. Using a series of controlled mixture experiments, we demonstrate that the eDNA signal from production water shows a strong, seasonally consistent linear relationship with herring weight fractions, however, the relationship is influenced by the molecular method used (qPCR or metabarcoding). In four large sprat landings analyzed, despite examples of remarkable consistency between eDNA and visual reporting, estimates of herring bycatch biomass varied between the methods applied, with the eDNA-based estimates having the highest precision for all landings analyzed. The eDNA-based bycatch assessment method has the potential to improve the quality and cost effectiveness of bycatch assessment in large pelagic fisheries catches and in the long run lead to more sustainable management of pelagic fish as a precious marine resource.

Genetic assignment predicts depth of benthic settlement for 0-group Atlantic cod.

Atlantic cod is a keystone species that remains among the most economically important demersal fish in the North Atlantic. Throughout its distribution range, Atlantic cod is composed of populations with varying environmental preferences and migratory propensities. This life-history variation is likely to have contributed to the niche width and large population sizes of Atlantic cod, and its relative resilience to environmental change and exploitation. The Icelandic cod stock is currently managed as a single unit, but early research indicates population variation by depth and temperature and distinct offshore and inshore spawning components. Pelagic 0-group juveniles from different spawning grounds coexist in nursery areas around Iceland, but their genetic composition or habitat partitioning had not been examined post benthic settlement. In the current study we examine the genetic composition of Atlantic cod juvenile aggregations at nearshore nursery grounds in NW-Iceland and report distinct segregation by the depth of offshore and inshore juvenile cod. The physiological mechanism of this segregation is not known, but the pattern demonstrates the need to consider population structure at nursery grounds in the application of marine spatial planning and other area-based conservation tools.

First tagging data on large Atlantic bluefin tuna returning to Nordic waters suggest repeated behaviour and skipped spawning.

Atlantic bluefin tuna (Thunnus thynnus; ABFT) is one of the most iconic fish species in the world. Recently, after being very rare for more than half a century, large bluefin tunas have returned to Nordic waters in late summer and autumn, marking the return of the largest predatory fish in Nordic waters. By tagging 18 bluefin tunas with electronic tags (pop-up satellite archival tags), we show that bluefin tuna observed in Nordic waters undertake different migration routes, with individuals migrating into the western Atlantic Ocean, while others stay exclusively in the eastern Atlantic and enter the Mediterranean Sea to spawn. We additionally present evidence of possible skipped spawning inferred from behavioural analyses. In Nordic waters, ABFT are primarily using the upper water column, likely reflecting feeding activity. The results support the hypothesis that ABFT migrating to Nordic waters return to the same general feeding area within the region on an annual basis. These observations may have important implications for management because (1) tunas that come into Nordic waters might represent only a few year classes (as evidenced by a narrow size range), and thus may be particularly vulnerable to area-specific exploitation, and (2) challenge the assumption of consecutive spawning in adult Atlantic bluefin tuna, as used in current stock assessment models. Without careful management and limited exploitation of this part of the ABFT population, the species' return to Nordic waters could be short-lived.

Fishing for DNA? Designing baits for population genetics in target enrichment experiments: Guidelines, considerations and the new tool supeRbaits.

Targeted sequencing is an increasingly popular next-generation sequencing (NGS) approach for studying populations that involves focusing sequencing efforts on specific parts of the genome of a species of interest. Methodologies and tools for designing targeted baits are scarce but in high demand. Here, we present specific guidelines and considerations for designing capture sequencing experiments for population genetics for both neutral genomic regions and regions subject to selection. We describe the bait design process for three diverse fish species: Atlantic salmon, Atlantic cod and tiger shark, which was carried out in our research group, and provide an evaluation of the performance of our approach across both historical and modern samples. The workflow used for designing these three bait sets has been implemented in the R-package supeRbaits, which encompasses our considerations and guidelines for bait design for the benefit of researchers and practitioners. The supeRbaits R-package is user-friendly and versatile. It is written in C++ and implemented in R. supeRbaits and its manual are available from Github: https://github.com/BelenJM/supeRbaits.

Atlantic herring (Clupea harengus) population structure in the Northeast Atlantic Ocean.

The Atlantic herring Clupea harengus L has a vast geographical distribution and a complex population structure with a few very large migratory units and many small local populations. Each population has its own spawning ground and/or time, thereby maintaining their genetic integrity. Several herring populations migrate between common feeding grounds and over-wintering areas resulting in frequent mixing of populations. Thus, many herring fisheries are based on mixed populations of different demographic status. In order to avoid over-exploitation of weak populations and to conserve biodiversity, understanding the population structure and population mixing is important for maintaining biologically sustainable herring fisheries. The aim of this study was to investigate the genetic population structure of herring in the Faroese and surrounding waters, and to develop genetic markers for distinguishing between four herring management units (often called stocks), namely the Norwegian spring-spawning herring (NSSH), Icelandic summer-spawning herring (ISSH), North Sea autumn-spawning herring (NSAH), and Faroese autumn-spawning herring (FASH). Herring from the four stocks were sequenced at low coverage, and single nucleotide polymorphisms (SNPs) were called and used for population structure analysis and individual assignment. An ancestry-informative SNP panel with 118 SNPs was developed and tested on 240 individuals. The results showed that all four stocks appeared to be genetically differentiated populations, but at lower levels of differentiation between FASH and ISSH than the other two populations. Overall assignment rate with the SNP panel was 80.7%, and agreement between the genetic and traditional visual assignment was 75.5%. The NSAH and NSSH samples had the highest assignment rate (100% and 98.3%, respectively) and highest agreement between traditional and genetic assignment methods (96.6% and 94.9%, respectively). The FASH and ISSH samples had substantially lower assignment rates (72.9% and 51.7%, respectively) and agreement between traditional and genetic methods (39.5% and 48.4%, respectively).

Replicated anthropogenic hybridisations reveal parallel patterns of admixture in marine mussels.

Human-mediated transport creates secondary contacts between genetically differentiated lineages, bringing new opportunities for gene exchange. When similar introductions occur in different places, they provide informally replicated experiments for studying hybridisation. We here examined 4,279 Mytilus mussels, sampled in Europe and genotyped with 77 ancestry-informative markers. We identified a type of introduced mussels, called "dock mussels," associated with port habitats and displaying a particular genetic signal of admixture between M. edulis and the Mediterranean lineage of M. galloprovincialis. These mussels exhibit similarities in their ancestry compositions, regardless of the local native genetic backgrounds and the distance separating colonised ports. We observed fine-scale genetic shifts at the port entrance, at scales below natural dispersal distance. Such sharp clines do not fit with migration-selection tension zone models, and instead suggest habitat choice and early-stage adaptation to the port environment, possibly coupled with connectivity barriers. Variations in the spread and admixture patterns of dock mussels seem to be influenced by the local native genetic backgrounds encountered. We next examined departures from the average admixture rate at different loci, and compared human-mediated admixture events, to naturally admixed populations and experimental crosses. When the same M. galloprovincialis background was involved, positive correlations in the departures of loci across locations were found; but when different backgrounds were involved, no or negative correlations were observed. While some observed positive correlations might be best explained by a shared history and saltatory colonisation, others are likely produced by parallel selective events. Altogether, genome-wide effect of admixture seems repeatable and more dependent on genetic background than environmental context. Our results pave the way towards further genomic analyses of admixture, and monitoring of the spread of dock mussels both at large and at fine spacial scales.

Northern European Salmo trutta (L.) populations are genetically divergent across geographical regions and environmental gradients.

The salmonid fish Brown trout is iconic as a model for the application of conservation genetics to understand and manage local interspecific variation. However, there is still scant information about relationships between local and large-scale population structure, and to what extent geographical and environmental variables are associated with barriers to gene flow. We used information from 3,782 mapped SNPs developed for the present study and conducted outlier tests and gene-environment association (GEA) analyses in order to examine drivers of population structure. Analyses comprised >2,600 fish from 72 riverine populations spanning a central part of the species' distribution in northern Europe. We report hitherto unidentified genetic breaks in population structure, indicating strong barriers to gene flow. GEA loci were widely spread across genomic regions and showed correlations with climatic, abiotic and geographical parameters. In some cases, individual loci showed consistent GEA across the geographical regions Britain, Europe and Scandinavia. In other cases, correlations were observed only within a sub-set of regions, suggesting that locus-specific variation was associated with local processes. A paired-population sampling design allowed us to evaluate sampling effects on detection of outlier loci and GEA. Two widely applied methods for outlier detection (pcadapt and bayescan) showed low overlap in loci identified as statistical outliers across sub-sets of data. Two GEA analytical approaches (LFMM and RDA) showed good correspondence concerning loci associated with specific variables, but LFMM identified five times more statistically significant associations than RDA. Our results emphasize the importance of carefully considering the statistical methods applied for the hypotheses being tested in outlier analysis. Sampling design may have lower impact on results if the objective is to identify GEA loci and their population distribution. Our study provides new insights into trout populations, and results have direct management implications in serving as a tool for identification of conservation units.

Genetic analyses reveal complex dynamics within a marine fish management area.

Genetic data have great potential for improving fisheries management by identifying the fundamental management units-that is, the biological populations-and their mixing. However, so far, the number of practical cases of marine fisheries management using genetics has been limited. Here, we used Atlantic cod in the Baltic Sea to demonstrate the applicability of genetics to a complex management scenario involving mixing of two genetically divergent populations. Specifically, we addressed several assumptions used in the current assessment of the two populations. Through analysis of 483 single nucleotide polymorphisms (SNPs) distributed across the Atlantic cod genome, we confirmed that a model of mechanical mixing, rather than hybridization and introgression, best explained the pattern of genetic differentiation. Thus, the fishery is best monitored as a mixed-stock fishery. Next, we developed a targeted panel of 39 SNPs with high statistical power for identifying population of origin and analyzed more than 2,000 tissue samples collected between 2011 and 2015 as well as 260 otoliths collected in 2003/2004. These data provided high spatial resolution and allowed us to investigate geographical trends in mixing, to compare patterns for different life stages and to investigate temporal trends in mixing. We found similar geographical trends for the two time points represented by tissue and otolith samples and that a recently implemented geographical management separation of the two populations provided a relatively close match to their distributions. In contrast to the current assumption, we found that patterns of mixing differed between juveniles and adults, a signal likely linked to the different reproductive dynamics of the two populations. Collectively, our data confirm that genetics is an operational tool for complex fisheries management applications. We recommend focussing on developing population assessment models and fisheries management frameworks to capitalize fully on the additional information offered by genetically assisted fisheries monitoring.

Parallel evolution and adaptation to environmental factors in a marine flatfish: Implications for fisheries and aquaculture management of the turbot (Scophthalmus maximus).

Unraveling adaptive genetic variation represents, in addition to the estimate of population demographic parameters, a cornerstone for the management of aquatic natural living resources, which, in turn, represent the raw material for breeding programs. The turbot (Scophthalmus maximus) is a marine flatfish of high commercial value living on the European continental shelf. While wild populations are declining, aquaculture is flourishing in southern Europe. We evaluated the genetic structure of turbot throughout its natural distribution range (672 individuals; 20 populations) by analyzing allele frequency data from 755 single nucleotide polymorphism discovered and genotyped by double-digest RAD sequencing. The species was structured into four main regions: Baltic Sea, Atlantic Ocean, Adriatic Sea, and Black Sea, with subtle differentiation apparent at the distribution margins of the Atlantic region. Genetic diversity and effective population size estimates were highest in the Atlantic populations, the area of greatest occurrence, while turbot from other regions showed lower levels, reflecting geographical isolation and reduced abundance. Divergent selection was detected within and between the Atlantic Ocean and Baltic Sea regions, and also when comparing these two regions with the Black Sea. Evidence of parallel evolution was detected between the two low salinity regions, the Baltic and Black seas. Correlation between genetic and environmental variation indicated that temperature and salinity were probably the main environmental drivers of selection. Mining around the four genomic regions consistently inferred to be under selection identified candidate genes related to osmoregulation, growth, and resistance to diseases. The new insights are useful for the management of turbot fisheries and aquaculture by providing the baseline for evaluating the consequences of turbot releases from restocking and farming.

Genetic diversity and connectivity within Mytilus spp. in the subarctic and Arctic.

Climate changes in the Arctic are predicted to alter distributions of marine species. However, such changes are difficult to quantify because information on present species distribution and the genetic variation within species is lacking or poorly examined. Blue mussels, Mytilus spp., are ecosystem engineers in the coastal zone globally. To improve knowledge of distribution and genetic structure of the Mytilus edulis complex in the Arctic, we analyzed 81 SNPs in 534 Mytilus spp. individuals sampled at 13 sites to provide baseline data for distribution and genetic variation of Mytilus mussels in the European Arctic. Mytilus edulis was the most abundant species found with a clear genetic split between populations in Greenland and the Eastern Atlantic. Surprisingly, analyses revealed the presence of Mytilus trossulus in high Arctic NW Greenland (77°N) and Mytilus galloprovincialis or their hybrids in SW Greenland, Svalbard, and the Pechora Sea. Furthermore, a high degree of hybridization and introgression between species was observed. Our study highlights the importance of distinguishing between congener species, which can display local adaptation and suggests that information on dispersal routes and barriers is essential for accurate predictions of regional susceptibility to range expansions or invasions of boreal species in the Arctic.

Historical DNA documents long-distance natal homing in marine fish.

The occurrence of natal homing in marine fish remains a fundamental question in fish ecology as its unequivocal demonstration requires tracking of individuals from fertilization to reproduction. Here, we provide evidence of long-distance natal homing (>1000 km) over more than 60 years in Atlantic cod (Gadus morhua), through genetic analysis of archived samples from marked and recaptured individuals. Using a high differentiation single-nucleotide polymorphism assay, we demonstrate that the vast majority of cod tagged in West Greenland and recaptured on Icelandic spawning grounds belonged to the Iceland offshore population, strongly supporting a hypothesis of homing. The high degree of natal fidelity observed provides the evolutionary settings for development of locally adapted populations in marine fish and emphasize the need to consider portfolio effects in marine fisheries management strategies.

Archived DNA reveals fisheries and climate induced collapse of a major fishery.

Fishing and climate change impact the demography of marine fishes, but it is generally ignored that many species are made up of genetically distinct locally adapted populations that may show idiosyncratic responses to environmental and anthropogenic pressures. Here, we track 80 years of Atlantic cod (Gadus morhua) population dynamics in West Greenland using DNA from archived otoliths in combination with fish population and niche based modeling. We document how the interacting effects of climate change and high fishing pressure lead to dramatic spatiotemporal changes in the proportions and abundance of different genetic populations, and eventually drove the cod fishery to a collapse in the early 1970s. Our results highlight the relevance of fisheries management at the level of genetic populations under future scenarios of climate change.

Spatiotemporal SNP analysis reveals pronounced biocomplexity at the northern range margin of Atlantic cod Gadus morhua.

Accurate prediction of species distribution shifts in the face of climate change requires a sound understanding of population diversity and local adaptations. Previous modeling has suggested that global warming will lead to increased abundance of Atlantic cod (Gadus morhua) in the ocean around Greenland, but the dynamics of earlier abundance fluctuations are not well understood. We applied a retrospective spatiotemporal population genomics approach to examine the temporal stability of cod population structure in this region and to search for signatures of divergent selection over a 78-year period spanning major demographic changes. Analyzing >900 gene-associated single nucleotide polymorphisms in 847 individuals, we identified four genetically distinct groups that exhibited varying spatial distributions with considerable overlap and mixture. The genetic composition had remained stable over decades at some spawning grounds, whereas complete population replacement was evident at others. Observations of elevated differentiation in certain genomic regions are consistent with adaptive divergence between the groups, indicating that they may respond differently to environmental variation. Significantly increased temporal changes at a subset of loci also suggest that adaptation may be ongoing. These findings illustrate the power of spatiotemporal population genomics for revealing biocomplexity in both space and time and for informing future fisheries management and conservation efforts.

Identification of single nucleotide polymorphisms in candidate genes for growth and reproduction in a nonmodel organism; the Atlantic cod, Gadus morhua.

Recent technological developments have facilitated intensified searches for genetic markers under selection in nonmodel species. Here, we present an approach for the identification of candidate gene variation in nonmodel organisms. We report on the characterization of 82 single nucleotide polymorphisms (SNPs) and on the development of a specific genotyping assay for 30 SNPs in 18 candidate genes for growth and reproduction in Atlantic cod (Gadus morhua). These markers can be used for scanning natural populations for signatures of selection in both contemporary and archived historical samples, for example in retrospective studies assessing the effects of environmental changes, such as increasing temperatures, and selection imposed by high fishing pressure. Furthermore, these gene markers may be of interest to aquaculture, serving as a starting point for linking phenotypic traits important for productivity with genotypes and potentially be of use for marker-assisted selection in the future. This study demonstrates that the candidate gene approach is a valuable and targeted complement to the more random approach for discovering genetic variation in the genome and transcriptome applied through high throughput methods in nonmodel species.

Evolutionary mechanisms shaping the genetic population structure of marine fishes; lessons from the European flounder (Platichthys flesus L.).

A number of evolutionary mechanisms have been suggested for generating low but significant genetic structuring among marine fish populations. We used nine microsatellite loci and recently developed methods in landscape genetics and coalescence-based estimation of historical gene flow and effective population sizes to assess temporal and spatial dynamics of the population structure in European flounder (Platichthys flesus L.). We collected 1062 flounders from 13 localities in the northeast Atlantic and Baltic Seas and found temporally stable and highly significant genetic differentiation among samples covering a large part of the species' range (global F(ST) = 0.024, P < 0.0001). In addition to historical processes, a number of contemporary acting evolutionary mechanisms were associated with genetic structuring. Physical forces, such as oceanographic and bathymetric barriers, were most likely related with the extreme isolation of the island population at the Faroe Islands. A sharp genetic break was associated with a change in life history from pelagic to benthic spawners in the Baltic Sea. Partial Mantel tests showed that geographical distance per se was not related with genetic structuring among Atlantic and western Baltic Sea samples. Alternative factors, such as dispersal potential and/or environmental gradients, could be important for generating genetic divergence in this region. The results show that the magnitude and scale of structuring generated by a specific mechanism depend critically on its interplay with other evolutionary mechanisms, highlighting the importance of investigating species with wide geographical and ecological distributions to increase our understanding of evolution in the marine environment.